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Abstract: We present results from quenched lattice QCD for the form factors for the 
decay B plv. The calculations are performed using a nonperturbatively improved action 
and operators at two values of the lattice spacing. The bottom quark mass is reached by 
extrapolation from simulations performed with heavy quark masses around the charm 
mass. Our primary result is for the partially integrated decay rate Fpi over the range 
12.7 GeV^ < g2 < 18.2 GeV^: 

rpi = 4.9l}^+^°xl0i2^-MKbp. 
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1. Introduction 

One of the primary goals of modern particle physics experiments is to determine the ele- 
ments of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. This matrix is unitary in the 
Standard Model (SM) and so-called unitarity triangles can be constructed using the or- 
thogonality of pairs of rows or columns. For the most common triangle, drawn in the 
complex p~rj plane of the Wolfenstein parameterisation of the CKM matrix , BaBar and 
Belle have produced the most accurate measurements of one of the angles (via sin 2/3) to 
date @, §]. To fully constrain the triangle one needs to know the length of the opposing 
side, governed by Vub- Any process in which a b quark transforms to a u quark can be 
used to measure V^b- Two possibilities are the exclusive decays B — > Trlv and B plv. 
The form factors which parameterise the hadronic amplitudes for these decays can be de- 
termined from first principles in lattice QCD. Combining these with measurements of the 
decay rate would allow the determination of V^b- In particular, combining a measured 
differential decay rate with the lattice values for the form factors at high would allow 
Vub to be extracted independently of any model for the behaviour of the form factors. 

The transition amplitude for semileptonic decays factorises into leptonic and hadronic 
parts when <C ?ti^, where q^ is the momentum transfer squared and my/ is the mass 
of the W boson. The hadronic matrix element contains the non-perturbative strong inter- 
action effects and is the largest source of uncertainty in theoretical determinations of the 
decay rates. 
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The hadronic matrix elements for B plv are parameterised by four form factors 

{p{k),i^\V,\B{p)) = .^Y^fLe .pPk-r^*' (1.1) 
ms + rup 

{p{k),v\A^\Bip)) = i{mB + mp)Ai{q^)gp^r^r -^^^^(.P + k)^q.vr 

niB + nip 

^2nrv^ (1.2) 

where B is the initial state pseudoscalar meson with three-momentum p, p is the final state 
vector meson with three-momentum k and polarisation rj, and q = p — k. The form factor 
A can be written as 

Aiq^) = Aoiq'') - Asiq^) (1.3) 

where 

= !I^^A,iq^) - ^^^A.iq^) (1.4) 
znip ZTJip 

with Ao{0) = AsiO). 

In the limit of zero lepton mass the differential decay rate is given by 0] 



dq^ 1927r3mf 



"B 

where 




H\q^) = (i^^^>l,(g2) - (ml - - g^) (m^ + rUp) A,{q 



2mp^q'^ [ms + mp 
HHi') = {ruB + mp) A,{q^) ± ^"'f ^(g^). 
Gp is the Fermi constant and A is the kinematic factor, 

\{q^)={ml+ml-q^f -AmWp (1-8) 

can be extracted from the experimental decay rate once the theoretical rate has been 
determined. The form factors Hq and H± correspond to the contributions of longitudinally 
and transversely polarised light vector mesons respectively 

We present results for the form factors of the decays B — > plv determined in the 
quenched approximation to lattice QCD, at two values of the lattice spacing. We use the 
non-perturbatively (NP) improved Sheikholeslami-Wohlert (SW) action and improved 
operators. Lattice artefacts then appear formally at O(a^) (where a is the lattice spacing) 
rather than 0{a), though this does not guarantee that they are numerically small. 

Mass-dependent lattice artefacts depending on [amq)'^ in the NP improved renormali- 
sation scheme limit the heavy quark mass mq at which one can simulate, without resorting 
to an effective action such as the Fermilab formalism Q. We have simulated heavy quarks 
with several different masses around the charm quark mass, then used continuum Heavy 
Quark Symmetry (HQS) to motivate the form of our extrapolation to the h quark mass. 
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In this way amg and thus {aruQ)^ are less than one. We have already presented a deter- 
mination of the form factors for B — > irlv at our finest value of the lattice spacing j^. 

The rest of the paper is organised as follows. Section two describes the details of the 
calculation and how the form factors were extracted from the data and then interpolated 
and extrapolated to physical quark masses. Section three discusses the results and contains 
an analysis of systematic errors. 

The results reported here represent a new analysis of the data used for an earlier 
preliminary analysis in 

2. Details of the calculation 

2.1 Improved action and operators 

In the Wilson formulation of lattice QCD, the fermionic part of the action has lattice 
artefacts of 0(a), while the gauge action differs from the continuum Yang-Mills action by 
terms of 0{a^). To leading order in a the Symanzik improvement program involves adding 
the SW term to the fermionic Wilson action, 

SssN = 5'w - csw^-^^^'4^{x)ia^yF^y{x)'4){x) (2.1) 

X 

Full 0{a) improvement of on-shell matrix elements also requires that the currents are 
suitably improved. The improved vector current is 

= V^{x) + acyduT^uix) 
Alix) = A^ix) + acAd^Pix) (2.2) 

where 

Vf,{x) = i){x)-i^'4}{x) 
A^{x) = i}{x)jf,'j5'lp{x) 

P{x) = ■ip{x)j5ip{x) 
T^u{x) = 'ip{x)ia^utp{x) 

and is the symmetric lattice derivative. The current renormalisation is as follows (J = 

A,vy. 

= Zj{l + bjamg)J^ (2.3) 

where Zj is calculated in a mass-independent renormalisation scheme. 
The bare quark mass, arug, is 

amg = U-- (2.4) 

where n is the hopping parameter. For non-degenerate currents, an effective quark mass is 
used in the definition of the renormalised current, corresponding to 

KcS 2 V^l K2J 
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In this renormahsation scheme, the improved quark mass, used in the chiral extrapolations, 
is defined as 



VHa 



mail + b^anin 



(2.6) 



2.2 Simulation details 



Gauge field configurations were generated using a combination of the over-relaxed |1C, 11| 
and the Cabibbo-Marinari |]12| algorithms with periodic boundary conditions at two values 



of the gauge coupling [5 = Q/qq. At each j3, heavy quark propagators were computed at 
four values of the hopping parameter, corresponding to quarks with masses in the region 
of the charm quark mass. For light quarks, three values of k were used for the light quark 
which occurs in the current, the active quark (A), and, owing to disk space constraints, 
only two values of k for the passive quark (P). All the light quarks had masses around that 
of the strange quark. Table |^ lists the input parameters. Kcrit, the value of the hopping 



/3 = 6.2 



/? = 6.0 



Volume 

csw 

-^configs 

(GeV) 
Heavy n 
Light K 



24^ X 48 
1.614 
216 
2.9i;j 

0.1200,0.1233,0.1266,0.1299 
0.1346,0.1351,0.1353 
0.13581+1 
0.135781^ 
0.134951^ 



16^ X 48 
1.769 
305 
2.12t} 

0.1123,0.1173,0.1223,0.1273 
0.13344,0.13417,0.13455 

0.13525li 

+1 



0.13520_^ 
0.13401+H 



Table 1: Input and derived parameters. The lattice spacing is set by tq. 



parameter which corresponds to zero quark mass. Kg and k„ are taken from |13] where 
the lattice spacing has been fixed from the Sommer scale tq |1^, 15|. Statistical errors are 



estimated using the bootstrap [16| with 1000 re-samplings. Unless otherwise specified, all 
errors listed in tables are statistical. 

Non-perturbative determinations of the improvement coefficients are available from 



two groups. The ALPHA collaboration have determined the value of csw llH using 
chiral symmetry and Ward identities in the Schrodinger Functional (SF) formalism. They 



have determined ca [18] and Za, Zy and by |19] in the same scheme and have a preliminary 

21]. Bhattacharya et al. |22, [2^, 24] have determined all the im- 



determination of cy 



provement coefficients needed to improve and renormalise quark bilinears, also using Ward 
identities, but on a periodic lattice with standard sources. They also use the ALPHA value 
of csw to improve the action. The values of Zy, Za and hy are in good agreement between 
the two groups. However, the ALPHA values for cy and ca are much larger than those 
of Bhattacharya et al. In a more detailed discussion of the improvement coefficients [25] 
we noted that the matrix elements determining leptonic decays are particularly sensitive 



to the value of the mixing c coefficients. As shown in equation (2^) it is the derivative 
of a current which is the improvement term. For the decay constants calculated at zero 
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three-momentum the derivative is the mass of the state which is 0{1). For the matrix 
elements determined here, the momentum component is a momentum difference and is 
0(1/10). Thus these matrix elements are relatively insensitive to the values of the mixing 
coefficients. We use the value of Bhattacharya et al. The value of bm is determined in 
lattice perturbation theory [26| with the boosted coupling = Qq/uq The mean link, 
uq, is taken from the plaquette expectation value, Uq = (ReTrC/p)/3. The values of the 
improvement coefficients are shown in table |2|. 

2.3 Correlation functions 



Details of the fitting procedure, and extrac- 
tion of the spectrum and amplitudes between 
the vacuum and meson states at zero three- 
momentum: (Mp I ^/;75'i/' 1 0) for example can be 
found in [|l^, . We used the fuzzing smear- 
ing function ||2^ for the light quark propaga- 
tors and a hydrogenic gauge invariant func- 
tion [29 1 for the heavy quarks. 

The light-light 2-point correlation func- 
tions at non-zero momentum were somewhat 
noisy. To determine these as accurately as 
possible we used a constrained dispersion re- 
lation fit. That is, we constrained the energy 





(3 = 6.2 


p = 6.0 


csw 


1.614 


1.769 


Za 


0.818(2)(5) 


0.807(2)(8) 


bA 


1.32(3)(4) 


1.28(3) (4) 


CA 


-0.032(3)(6) 


-0.037(4)(8) 


Zv 


0.7874(4) 


0.770(1) 


bv 


1.42(1)(1) 


1.52(1) 


cv 


-0.09(2)(1) 


-0.107(17)(4) 


bm 


-0.6517 


-0.6621 



Table 2: Improvement coefficients used in this 
work. 



of the non-zero momentum states to be 



(2.7) 



Furthermore, we constrained the amplitudes by fitting the correlator which was fuzzed at 
source and local at sink (FL) to the following form 



C{p,t) 



Z^{\p?)Z' 
2E 



-Et 



(2.8) 



where T is the time extent of the lattice and E is constrained from equation (|2.7| ). The 
local Z^, by Lorentz invariance, is not a function of momentum and so both and M 
are determined at zero momentum. 

We determined the amplitudes and energies of the heavy-light mesons at all momenta 
using a free single exponential fit. We then compared to the dispersion relation as a check 
of lattice artefacts. 

We have extracted the matrix element for the semileptonic decays from three-point 
and two-point correlation functions. The three-point function for the transition from state 
A to state B is given by, 

^e-*(^-^+(^-P^-'^)(Ol T{f)B(x)J^(2/)!^|4(0)} |0) 

x,y 



(2. 



where 0^ is the operator which creates state A, fi^ is the operator which destroys state B 
and is the weak current. We have computed the correlation function using the "standard 
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source" method po| , 31 1. We chose the extension timeshce to be, tx = 28 rather than in 
the centre of the lattice at t = 24. This enabled us to look for different systematic effects 
by comparing both sides of the lattice, for example to look at different time orderings of 
operators and to check the quality of plateaux. However, we lost the ability to average 
results from the front and back sides of the lattice and so our statistical errors are larger. 
When the operators are well separated, that is when ty and — ty are large, then 

C^^,ip,tx,k,ty) = ——ZA^^^Zl{p{k)\J>^mB{p)) (2.10) 

where Zi = (0|Pi|Pj(pi)). Both the Z's (for smeared operators) and the energies are 
functions of three- momentum squared, and can be extracted from the appropriate 

two-point functions, as described above. The remaining matrix element is the one in 
equation ([L1| ). 

The simulation used heavy meson spatial momentum of magnitude or 1, in lattice 
units of 27r/aL. The light meson was given spatial momentum of or 1 in lattice units. 
Six momentum channels with different values of were considered, shown in table |3[ 



channel 


P 


k 











(0,0,0) 


(0,0,0) 




1 





1 


(0,0,0) 


(1,0,0) 


+ perms 


5 


1 


2 


(1,0,0) 


(1,0,0) 




3 


2 


3 


(1,0,0) 


(0,0,0) 


+ perms 


4 


2 


4 


(1,0,0) 


(0,1,0) 


+ perms 


10 


3 


5 


(1,0,0) 


(-1,0,0) 




5 


1 



Table 3: Values of three-nioinentum in lattice units of and the number of distinct matrix 
elements for each momentum channel. "+ perms" implies permutations of k that have the same 



For each channel, different spatial and temporal components combined with the rel- 
evant three-momentum allow equations ( |1.1[) to be decomposed into distinct matrix el- 
ements. The three-point correlation functions were fitted simultaneously to the distinct 
matrix elements with the form factors as free parameters. We averaged over momenta 
and Lorentz channels that have the same matrix element. The number of distinct matrix 
elements for each channel and current is shown in Table ^. For channel 0, the coefficients 
vanish for all the form factors except Ai. 

Figure |l] shows examples of the averaged three-point correlation functions with the 
time dependence divided out. The data on the fore-side of the lattice has larger statistical 
errors. This is simply due the noise-to-signal ratio increasing at larger time separations. 
We observed clear plateaux, but note that the plateaux on the back-side of the lattice 
are better for extracting the signal. We chose to fit to the three-point function with the 
time-dependence cancelled using fitted two point functions rather than use the ratio of 
three-point over two-point correlation functions as the current operator has less overlap 
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with excited states. This is supported by the data, as we see longer plateaux when using 
the fitted two-point parameters. 

It is clear that even on the shorter back- 
side of the lattice the separation of the op- 
erators is sufficient to see only the ground 
state. For the spatial current (open sym- 
bols) there is a significant disagreement be- 
tween each side of the lattice. This may 
be due to contamination from different time 
orderings of the operators wrapping around 
the lattice. In this case, one would expect 
the shorter back-side of the lattice to be less 
affected than the longer fore-side. 

In general, we fitted to the central five 
timeslices on the back side of the lattice. 
For several channels we have shifted this 
window to accommodate a better plateau. 
Some channels had much clearer plateaux 
on the fore-side of the lattice. In all cases 
the criteria for choosing the fit range was 
a plateau such that the form factors were 
stable against a change of one timeslice in 

the fit range and that had a reasonable reduced and Q-value. We used a block diag- 
onal covariance matrix to prevent different systematics between different currents being 
misinterpreted as correlation information. 

By measuring the quark mass dependence of the form factors, we can interpolate 
and/or extrapolate to the quark masses of the physical decays we are interested in. 

2.4 Light quark dependence 

In quenched lattice QCD calculations, light quarks are typically simulated around the 
strange quark mass and then extrapolated to light or zero quark mass. The reasons for 
this are mainly algorithmic and computational. In the first instance, the computation 
of light quark propagators suffers from critical slowing down. The number of iterations 
it takes to invert the fermion matrix varies inversely with the quark mass, making light 
quarks prohibitively expensive in computer time. Second, as the quarks become lighter 
the hadron states become larger, requiring a larger box size, again increasing the amount 
of computation required. Third, quenched QCD has the wrong chiral behaviour, but in 
practice it is the computational requirements which restrict the range of quark mass. 
The form factor F has both explicit, and implicit mass dependence. 



Figure 1: Cspt with measured two-point fac- 
tors divided off, for (3 = 6.2, light quarks, ka — 
Kp = 0.1346 and the heavy quark, k = 0.1200, 
kinematic channel 1. Different symbols show the 
distinct matrix elements for (V^|A|P). 



F = F {q'^{mA,mp),mA,mp) 



(2.11) 



where niA is the active quark mass and mp is the passive quark mass. For small changes 
in niA, rnp, and the variation in the form factor can be approximated by a first-order 
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Taylor expansion 



dF 



F{q^,mA,rnp) = a + mA-^ h mp-^ h q^-^ 

am A om p oq'^ 



dF ^dF 



(2.12) 



For a particular momentum channel the exact dependence of on meson mass is known. 
The first-order Taylor expansion of q^ gives 



q'^ = (3 + jMh + 6Ml 



(2.13) 



where Mh is the mass of the heavy pseudoscalar meson and Ml is the mass of the light 
state. To a good approximation Mh varies linearly with mp. The light vector meson mass 
depends linearly on the averaged light quark mass mes = (rriA + mp)/2. 

To control this extrapolation, we separated the chiral and momentum behaviour of 
the form factors by first interpolating in q'^, at fixed quark mass, to a common set of q'^ 
values over the range of quark masses. Then we extrapolated in quark mass at fixed g^, as 



suggested by To interpolate the form factors in q , we used pole models |3^, |35|] , 



F{q' 



m 



1 — q^ /m^p 



(2.14) 



This is illustrated in figure ^. As we have interpolated in q^ any model dependence will 
be small. We have checked for residual model dependence by varying the model. The 
figure shows the form factors and the pole mass fits. We note that before any heavy quark 
extrapolation we can cover the full range of physical q^ values. 



1.75 



1.25 



3 




0.75 



0.25 



Q I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

-0.1 -0.05 0.05 0.1 0.15 0.2 

2 2 

a q 



Figure 2: The form factors as a function of g^, at /3 = 6.2 for kh = 0.1200, k^i = Kp = 0.1346. 
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Once the dependence of the form factors had been determined, we chose a set of 
common (f' values to which we interpolated at each value of the quark mass. We then 
extrapolated the form factors at fixed in two dimensions according to 

F{mcs,mp) = a + /3mp + 7mefr (2-15) 

g2=const 

2.5 Heavy quark extrapolation 

Near zero-recoil (^^ax) form factors scale in the following way |36], 



F ~ (2.16) 

where = 1 for F G {^,^0,^2}, = —1 for F = Ai. We introduce a new kinematic 
variable 

vk= ^ — (2.17 

where v is the four-velocity of the heavy meson. We then have the following form for the 
extrapolation at fixed v ■ k. 

= C {Mh, Mb) F{v ■ m-""^ = c (1 + ^ + ^ + • • •) (2.18) 



The coefficient C is the logarithmic matching factor |37| 



C,M„,M,)H«)^'" (2.19) 



MS 



where /5o = H in quenched QCD and is the one- loop running coupling with 
295 MeV. 

With the fixed-g^ extrapolation, it is easy to implement the heavy quark extrapolation. 
We used a different set of q^ at each value of the heavy quark mass so that the heavy quark 
extrapolation could be done at fixed v ■ k. 

Figure ^ shows the extrapolation of the form factor Ai. The solid line shows the 
quadratic fit to all four data points at /3 = 6.2, and the dot-dashed line shows the linear 
fit to the heaviest three. At the B meson mass the two lines are still close together, 
suggesting that the extrapolation is reasonable. Comparing the circles (/? = 6.2) with the 
squares (/? = 6.0) there is little difference, suggesting that 0{{amf') lattice artefacts are 
small. In ||2^ we fitted simultaneously to the data for the decay constants at both lattice 
spacings, allowing for lattice artefacts, and found the effect to be small. The method of 
extrapolation is the same as used here. However, we do not attempt a simultaneous fit 
since the statistical precision of the data is not sufficiently good even for Ai . 

3. Results 

In principle one could compute the form factors for any value of q'^ from lattice QCD. 
However, states with high spatial momentum are very noisy and thus difficult to measure 



-9- 




on the lattice. Thus we are restricted to the high end of the range. In addition, the 
procedures we have introduced to control the extrapolations, separating the from the 
quark mass dependence, have further restricted the range of away from g^^x' the 
range 

12.7 GeV^ <q^ < 18.2 GeV^ (3.1) 

Moreover, the relatively small number of momentum channels for which the form factors are 
extracted, six for Ai, five for Aq and A2, and four for V ^ coupled with the interpolation at 
fixed q^ imply by naive counting of degrees of freedom that we have only four independent 
data for Ai and worse, two independent data for V . Fitting the functional form of the (f' 
dependence of the form factors is thus rather hard. However, we are free to evaluate the 
form factors, and thus the differential decay rate, at any value of q^ we choose without 
introducing any extra model dependence as long as it is in the range of allowed q^. In 
particular we can determine a partially integrated decay rate over this range. 

Figure Q shows the four form factors on both lattices. In this case we have chosen nine 
values of (f' . The form factor Ai which dominates at (/max is well determined and is in good 
agreement for both lattice spacings. The subleading form factors have a much noisier signal, 
especially for the coarser lattice. This made the extrapolations very difficult to control. 
For the coarse lattice only we introduced two additional model dependent constraints on 
the data during the fixed quark mass q^ interpolation. First, according to pole models the 
pole mass for Ai and A2 should be the same, and we enforce this constraint. As Ai is 
much better determined than A2, Ai remains unchanged. Second, we find a dipole rather 
than a pole fits the data better for V . Both these constraints affect the value of the form 
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factors very little for the fine lattice and are not required to control the extrapolations. As 
the finer lattice forms our result and the coarse lattice is used to estimate lattice artifacts, 
the resulting model dependence enters only in our estimate of systematic error. 

(GeV)^ (GeV)^ 



10 12 14 16 18 20 12 14 16 18 20„ 




Figure 4: The form factors on both lattices. The vertical scale is different for each form factor. 

In Table ^ we list the values of the form factors on both lattices at the highest and 
lowest values of q'^ only. Using equation ( |1.5D we can similarly determine the diff'erential 
decay rate over the same range of momentum. Shown in figure ^ is the differential decay 
rate on both lattices. To reiterate, we can evaluate the differential decay rate for any 
in the allowed range but only two data points will be independent. For this reason, 
our main result is for the partially integrated decay rate over the allowed range given in 
equation (|3.l|) 

712.7 GeV2 5r 

where the first error is statistical and the second is systematic. Estimates of systematic 
error are discussed in the next section. 

The CLEO collaboration [^] have measured the rates for the decays B ttIu and 
B — > pliy. They determine the partial decay rate for three ranges of g^: < q^ < SGeV'^, 
8 < g2 < l6GeV^ and q'^ > IGGeV^. They compare the measured partially integrated 
decay rate with Light Cone Sum Rules (LCSR) results for q"^ < 16 GeV^ and to lattice QCD 
data for > 16 GeV'^. We cannot bin our data in the same way as we have not determined 
the form factors at g^ax) cannot directly compare. The BaBar collaboration [^] have 

measured the total decay rate, but only quote for the lifetime and V^h, so again we cannot 
compare our results to the experimental data at the same values of q"^. 
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Figure 5: The differential decay rate on both lattices. The vertical line to the left is the charm 
end point, and to the right is 



1^ 

imax 



3.1 Systematic uncertainties 

This calculation is done in the quenched ap- 
proximation. There is therefore an uncontrolled 
systematic error. For most dimensionful quanti- 
ties such as the hadron spectrum this is of the 
order of 10%. The quenching error is manifest in 
the choice of quantity used to set the scale: in 
the quenched approximation different quantities 
such as the nucleon mass, or the Sommer scale, 
give different answers for the lattice spacing. The 
form factors are dimensionless and so are not di- 
rectly affected by the scale-setting. However there 
is an indirect dependence. We vary the choice of 
quantity used to set the scale. This in itself is 
not an estimate of the quenching systematic on 
the form factors, only of their implicit scale de- 
pendence. We use the Sommer scale, rg to set the 
scale. At each lattice spacing the value of a/rg is 
unambiguously defined although the experimental Table 4: The form factors on both lat- 
value is not known. Sommer originally advocated tices 
ro = 0.5 fm but determinations of the lattice spac- 



(GeV)2 


/J = 6.2 /3 = 6.0 




Ax 


18.2 
12.7 


0.60l^ 0.601 9 
0.45l^ 0.40l}^ 






18.2 
12.7 


0.881^°° 1.231,^0^3 
0.721 1 0.821 






18.2 
12.7 


2.28l^? 0.521216 
0.931}? 0.801 3I 




V 


18.2 
12.7 


9 rio+72 q 1 7+247 
^•^^-65 "^•-'^'-226 

0.741^0 0.461 1 




#T X IO12 s-^ 


18.2 
12.7 


1.4+ I l.Q+^l 
l.lt^l 0.6l \ 
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ing using the kaon decay constant, the nucleon mass or the rho mass correspond to rg values 
in the range 0.5-0.55 fm. We take ro = 0.5 fm as our central value and use rg = 0.55 fm to 
estimate the systematic error from quenching scale ambiguity. 

We simulate heavy quarks with several quark masses around the charm quark mass 
and then use continuum heavy quark symmetry to extrapolate the form factors to the 
bottom quark mass. We estimate the systematic uncertainty coming from the quadratic 
extrapolation in equation ( 2.1^ ) by performing a linear extrapolation to the heaviest three 



quark masses. This is shown as the dot-dashed line in figure ^ For the form factor Ai on 
the finer lattice spacing at high this effect is smaller than the statistical uncertainty. For 
the other form factors the effect is larger, so we include this in the systematic uncertainty 
of the partially integrated decay rate. 



In 1 25] we made several estimates of systematic uncertainty arising from the heavy- 
quark formalism and renormalisation method. We argue here that the same effects will be 
no larger for the form factors than for the heavy-light decay constants. We simultaneously 
fitted the heavy-light decay constant data at both lattice spacings and allowed for mass 
dependent lattice artefacts. This produced a 5% effect in the value of decay constant. We 
do not attempt this procedure here as our statistical errors are too large, but have no reason 



to believe it would be a greater effect for the form factors. As suggested by Bernard |40| 
the normalisation of the currents can be altered so that while formally the same at 0{a) 
it looks much more like the mass dependent normalisation of the Fermilab formalism [Q]. 



^ — > -0' = + fjam tp (3.3) 

where depends on the improvement coefficients Zj, cj and bj and the improved currents 
defined in equation (p.2|). Combining this alternative normalisation with the kinetic mass 
M2 defined as 

1 d'^E 



M2 dp 



(3.4) 

p=0 



we found that this produced an effect on the pseudoscalar decay constant smaller than the 
difference between the linear and quadratic extrapolations at the coarser lattice spacing. 
As we described in section |2.1| , the effect of the improvement for the form factors is much 
smaller than for the decay constants, so the overall effect of changing the normalisation 
and hadron mass cannot be greater than for the decay constants. Finally we varied the 



values of the improvement coefficients by the errors quoted by Bhattacharya et al. |22] and 
found the overall effect to be 1% in the axial and 3% in the vector currents. Again, as 
the effect of improvement is smaller in the form factors the overall effect cannot be greater 
than this. 

The main systematic uncertainties are tabulated in table |5|. To obtain the systematic 
error on the partially integrated decay rate we have combined the systematic errors in 
quadrature. 
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lattice spacing -20% -45% +14% 

quenched scale ambiguity -20% -27% -14% 

heavy quark extrapolation —16% —9% —14% 

pole model dependence —4% —5% +0% 

Table 5: Systematic uncertainties in the decay rate 
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